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ABSTRACT 

We use the Stripe 82 proper motion catalogue of |Bramich et al. ( 20081) to study the kinematics of 
Galactic disk stars in the solar neighborhood. We select samples of dwarf stars with reliable spectra 
and proper motions. They have cylindrical polar radius between 7 < R < 9 kpc, heights from the 
Galactic plane satisfying \z\ < 2 kpc and span a range of metallicities —1.5 < [Fe/H] < 0. We develop a 
method for calculating and correcting for the halo contamination in our sample using the distribution of 
rotational velocities. Two Gaussians representing disk and halo populations are used to fit the radial 
(vr) and vertical (v z ) velocity distributions via maximum likelihood methods. For the azimuthal 
velocities (v^) the same technique is used, except that a skewed non-Gaussian functional form now 
represents the disk velocity distribution. This enables us to compute the dispersions ctr, <t z , cr^ and 
cross-terms (the tilt ctr z and the vertex deviation ctr^) of the velocity ellipsoid as a function of height 
and metallicity. We also investigate the rotation lag of the disk, finding that the more metal-poor stars 
rotate significantly slower than the metal-rich stars. These samples provide important constraints on 
heating mechanisms in the Galactic disk and can be used for a variety of applications. We present one 
such application, employing the Jeans equations to provide a simple model of the potential close to 
the disk. Our model is in excellent agreement with others in the literature and provides an indication 
the disk, rather than the halo, dominates the circular speed at the solar neighborhood. We obtain 
a surface mass density within 1.1 kpc of around 66 M Q pc -2 and estimate a local halo density of 
0.015 M pc- 3 = 0.57 GeV cm- 3 . 

Subject headings: Galaxy: disk — Galaxy: evolution — Galaxy: kinematics and dynamics — solar 
neighborhood 



1. INTRODUCTION 

In the simplest model of a spiral galaxy, stars in 
the disk librate about circular orbits in t he equato 
rial plane ( "th e epicyclic approximation" , Binney & 



Tremaine||2008[ ). This picture is complicated by astro 
physical processes: such as bar instabilities, satellite ac 
cretion, scattering by spiral arms or collisions with molec 
ular clouds. Surveys which gather kinematic data allow 
astronomers to analyze the properties of large samples 
of stars and hence probe the signatures left by heating 
processes in their net motions and dispersions. For ex- 
ample, the Hipparcos satellite provided precise parallaxes 
and proper motions for nearby stars, which were used to 
provide the velocity dispersion as a functio n of stellar 
type in the im mediate solar neighborhood (Dehnen & 
Binney||1998a[ ). 



nant trapping, which converts planar orbits into inclined 
ones ( Sridhar fc Touma|| 19961), o r the dissipational col- 
lapse of gas ( Brook et al.| |2004), or the accretion of a 
satellite galaxy, in which case the thick disk may form 
from satellite debris, or a s a result of heating of the thin 
disk by the collision (eg. Villalobos & Helmi 2008). In 
recent years a lot of attention has tocussed on the torma- 
tion of thick disks through radial migration, where stars 
trapped in co-rotation resona nces of the spiral arms are 
transported through the disk (Scllw ood fe Binney||2002| 
Roskar et al.|2008 Schonrich fc Binney|2009p o r through 



similar~mecnamsms due to the Galactic bar (Minchev 



In their 1983 study of the density distribution of stars 
up to 4 kpc below the plane, Gilmore & Reid argue that 
the disk is in fact composed of two distinct components: 
a thin disk with approximate scale height ~ 300 pc, and 
an older, more metal-poor thick disk. Further studies 



& Famaey 2010). Detailed tests of the chemical prop- 
erties of trie Milky Way disk are also becoming possi- 
ble, for example using large samples of stars with alpha- 
element abundances to characterize the origins of the dif- 
ferent components (e.g. Bovy et al .|20 11 | [Lee et al.|2011 
Navarro et aL|201l[ [Ruchti et al.||20l0| ) 
Much can also be learnt about h< 



leatmg processes 

through studying the size and shape of the stellar veloc- 
ity ellipsoid, or equivalently the stellar velocity dispersion 
tensor. The ratio of radial to vertical velocity dispersion 
cr/cz, as well as the dependence of the dispersion on 
the age of a stellar population, can be compared with 
analytical and numerical predictions for various models 
fuTe^arose is'still opem Suggested scenario7incTude"reso- of disk , heat ing, to see which mech anisms dominate in 



(see |Reid fc Majewski|1993| for a review) have supported 
this claim, although th e thick disk scale height remains 



poorl y constrained (eg. Juric et al 2008 De Jong et al. 
2010 1 and the question of how the two-component struc 
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our galaxy ( |Fuchs fc Wieleri 19871. The ratio of tan- 
gential to radial velocity dispersion o^/ctr is predicted 
to be 1/ \/2 from epicylic theory for a galaxy with a flat 
rotation curve, though it has long been known to be less 
than this for the Milky Way ( |Kerr fc Lynden-Bell|[i~986 
Evans & Collett 1993). The covariances ctr z and or^ 
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are also of interest; non-zero ctr^ i mplies that the Galac- 
tic p otential is non-axisymmetric (Binney & Tremaine 
20081. The quantity ctr, z , or more specifically, its verti- 
cal gradient, is important for the calculation of Galactic 
paramete rs, including the local su rface mass density of 
the disk (jKuijkcn & Gilmorc 19 89]) and the asymmetric 
drift relation (Dehnen & Binney 1998a|. 

Large-scale photometric and kinematic surveys are 
making measurements of the components of the velocity 
dispersion tensor possible for populations in the Milky 
Way Galaxy. For example, u sing data from the R adial 
Velocity Experiment (RAVE) , [Siebert et al.| ( |2008[ ) mea- 
sured the magnitude and orientation of the velocity dis- 
persion tensor for red clump stars between 500 pc and 
1500 pc below the Galactic Plane. They reckoned the 
tilt to be 7.3 ± 1.8°, which is c onsistent with alignment 
in spherical polar coordinates. Siebert et al.| compare 
this value to computed inclinations for two mass mod- 
els of the Milky Way. The measurement is reproducible 
with a short scalelength of the stellar disk (« 2 kpc) if the 
dark halo is oblate or with a long scalelength (« 3 kpc) 
if the dark halo is spherical or prola te. A similar study 
has been published subsequently by Ca setti-Dinescu et] 



of dwarf stars. The procedure for splitting the sample 
into disk and halo components is described in Section 
[3] as well as a discussion of the effects of measurement 
errors. Our results are presented and discussed in Section 
® 

2. SAMPLE CONSTRUCTION 
We use data from the 7th public data release of the 



Sloan Digital Sky Survey (SDSS; Abazajian et al. 2009 



Yanny et al. 2009). We restrict attention to stars in the 



spectroscopic part of the survey, so as to have estimates 
of velocities, metallicities and surface gravities. The full 
details describing the analysis of these data, in particular 
the reliability of the de rived veloci t ies an d stellar param- 
eters, can be found in Lee et al. (2008) and references 
therein. 

Although SDSS data cover a huge part of the sky, 
mainly in the Northern Galactic cap, we also restrict 
ourselves to data in Stripe 82, so as to exploit the cat- 
alogue of high-precisio n photometry and proper-motions 



(Bramich et al. 2008). This stripe covers around 250 



al.| ( |2011 [), also using radial velocities from the RAVE 
survey, finding a tilt of 8.6 ± 1.8° for red clump stars at 
a distance of 700 to 2000 pc from the plane. 

The kinematics properties of the Galactic disk can 
also be used to probe directly its gravitationa l potential, 
as was carried out in the classical work of Kuij ken "fc] 
Gilmore ( 1991 ). Using a sample of around 500 K dwarfs, 
they were able to estimate that the surface mass density 
within 1.1 kpc was 71 ± 6 M pc -2 . Subsequent work 



has supported th is measurement, for example, Holmberg 
& Flynn (2004). An inventory of the surface density 
of baryonic material in the solar neighborhood indicates 
that around a third of this mass is likely non-baryonic 
(e-g " 



Binney & Evans 
be used to determine 



2001 1. Such techniques can also 
he local mass density and hence 



make predictions for the local dark matter de nsity (e.g. 
Holmberg fc Flynn|[2000]jGarbari et al.||2011[ ), which is 
a vital ingredie nt in predictions tor the"direct detection 
of dark matter ( Bertone, Hooper & Silk 2005 , and refer- 
ences therein). 

Clearly obtaining unbiased determinations of the kine- 
matic properties of the disk, including their gradients, is 
of great importance for understanding its structure and 
evolution. Here, we study the kinematics of disk stars us- 
ing data from the recently construc ted Sloan Digital Sky 
Survey (SDSS) Stripe 82 catalogue ( |Bramich et al.|2008| ). 
This covers an equatorial stripe of area 250 deg z , which 
has been repeatedly imaged by the SDSS since 1998 pri- 
marily with the aim of supernova discovery. The cata- 
logue contains almost 4 million "light-motion curves" and 
is complete down to a magnitude of r « 21.5. It reaches 
almost 2 mag fain ter than the SDSS/USNO-B catalogue 
( Munn et al. |2004 1 , making it the deepest large-area pho- 
tometric and astrometric catalogue available. We have 
already exploited the Stripe 82 catalogue to measure the 
components of th e velocity dispersion tensor for a sample 
of halo stars (e.g. Smith et aH2009a]b |. Here, we provide 
a complementary study for disk populations. 

In Section [2] we describe the selection cuts that we 
apply to the Stripe 82 catalogue to generate our sample deals with the analysis of photometric metal 



square degrees in region a = 20.7 hr to 3.3 hr and 
\S\ < 1.26°. This corresponds to the South Galactic cap, 
with I ~ 50 to 190 degrees and b ~ —25 to —60 degrees. 

From this data, we select stars with reliable spectra 
(i.e. labelled in the catalogue by the flag 'nnnn') and 
errors in log<? and [Fe/H] of less than 0.5 dex and in 
radial velocity of less than 20 kms -1 . We require proper 
motion errors of less than 4 mas yr -1 . We also place a cut 
in log g in order to remove giant stars from the sample, 
which we place at log g = 3.5. To investigate the effects 
of giant contamination, we repeat our analysis of the 
sample properties using the stricter selection criterion 
log g > 4. This stricter cut excludes only ~ 1 per cent 
of our sample. We checked that this made no qualitative 
differences to the trends observed. 

The errors quoted in the SDSS stellar parameter 
pipeline are internal errors. Comparison with obser- 
vations made using high-resolution spectroscopy reveals 
that there are also external errors of 2.4 kms -1 , 0.11 
dex and 0.21 dex present in the radial ve locities, metal- 
licities and surface gravities, respectively ( Allende Prieto 
et al.|[2008 |. We add the internal and external errors in 
quadrature to obtain the overall error in each measure- 
ment. 

We correct all magnitudes for the effec ts of extinction 
by int erstellar dust, using the maps of |Schlegel et al.| 
(1998). We also place a cut in color of 0.3 < {g-i) < 4.0. 
This is required so that we can obtain reliable photomet- 
ric distance s , whi ch we estimate using the relation of 
Ivezic et al. (2008 equations A1-A5 ) . We make o ne mi - 
nor modification to the relation o f|Ivezic et al.| ([2008), 
namely we do not apply the turn-off correction given by 
equation (A6) of their paper. Their correction is based 
on the metal-poor globular cluster M13, which may not 
be suitable for disk stars. In Appendix |A"] we use stellar 
models to show that a better choice may be to simply 
neglect this correction term. Note that although |Ivezic| 
et al. ( 2008 1 use photometric estimates for the metallic- 
ities in their paper, we are able to use the more reliable 
spectroscopic metallicities for our sample . It is worth 



pointing out that although the Ivezic et al. ( 2008 ) paper 



lcities, their 

parallax relation is derived from star clusters; as a con- 
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sequence our work will not be affected by any systematic 
offsets between spectroscopic and photometric metallic- 
ities (should any be present). The median error on our 
distances is 10.7 per cent. 

Once we have an estimate for the distance, we can 
determine the full three-dimensional positions and ve- 
locities. We express these in terms of Galactocentric 
cylindrical coordinates, (R, </>, z), where cj> is defined as 
increasing in the direction opposite to the solar rota- 
tion, and z is positive towards the North Galactic Pole 
(NGP). We calculate corresponding velocities (v r, v^, v, 



for each star, an d correct for the solar motion (Dehnen 
& Binney 1998a) and the motion of the local standard 
of rest (LSR) so that our stellar velocities are relative to 
the Galactic frame, ie. vlsr = — 220 km s -1 . 

The errors on the positions and velocities are calcu- 
lated using a Monte Carlo method. This is done by tak- 
ing the uncertainties on all of the observed quantities for 
a given star (proper motion, radial velocity, magnitudes 
and metallicity) and randomly sampling from these as- 
suming the errors are Gaussian; for each realisation we 
then calculate a distance (including uncertainties in the 
photometric parallax relation) and hence a velocity. We 
repeat this sampling 1000 times and from the resulting 
distributions of positions and velocities calculate our er- 
ror matrix for a given star. The median values of the 
errors 5vr, Sv^ and <5v z are 27.2, 27.3 and 16.9 kms^ 1 , 
respectively. As our line of sight is approximately aligned 
with the South Galactic Pole, v z is mostly determined 
by radial velocity measurements, whereas vr and v^ are 
mostly dependent on proper motions. The measurement 
errors in the proper motions are much larger than in 
the radial velocities, hence Svr and are larger than 
Sv z . We also calculate the correlations in the errors (e.g. 
5vrv z ), which are important for the covariances. 

3. ALGORITHMS FOR DISK KINEMATICS 

Here, we describe the techniques used to obtain the 
kinematic properties of the disk. We aim to study the 
properties of our sample as a function of height above 
or below the Galactic plane z. However, many of these 
properties depend strongly on metallicity, so to allow for 
the influence of vertical metallicity gradients in the disk, 
we split the sample into ranges in metallicity and study 
trends in each range separately. This is also crucial be- 
cause the SDSS spectroscopic selection function can in- 
troduce biases which are difficult to model. Therefore, 
by binning our data in metallicity ranges, we are only 
making the assumption that the metallicity distribution 
within a bin is representative of the true metallicity dis- 
tribution. Although this will not be exactly the case due 
to variations within a bin, both in terms of metallicity 
and age, it should remove the main source of systematic 
bias in the analysis. 

The range of our data in cylindrical polar radius R 
is too narrow for us to determine trends in the radial 
direction with any confidence. Clearly, the further we 
go from the plane, the larger the range in R which is 
covered, and so in order to keep our data to a specific 
range in R we restrict ourselves to 7 < R < 9. There- 
fore, we split the data into three ranges in metallicity 
(-1.5 < [Fe/H] < -0.8,-0.8 < [Fe/H] < -0.5 and 
—0.5 < [Fe/H]), and then for each metallicity bin we 
further divide the data into four ranges in z out to a 



maximum distance of 2 kpc. A total of 7280 stars match 
these criteria. The stars are equally divided between the 
four distance bins, resulting in in around 500 to 800 stars 
per bin (see Table [l}. 

3.1. The Stellar Halo Contribution 

In order to carry out an unbiased study of the kinemat- 
ics of the disk, we need a method to model the contami- 
nation from the halo stars in our sample. To do this, we 
make the assumption that all counter-rotating (v^ > 0) 
stars belong to the halo, i.e., the number of halo stars is 
simply twice the number of counter-rotating stars. This 
provides an estimate for the level of halo contamination 
in each bin, which allows us to make corrections when 
calculating our kinematic properties. Some auth ors have 
argued tha t the local stellar halo is rotating (e.g. Carollo 
et al. | 200"7 ), although this is in conflict with later studies 
e.g. [ivezTc et al." 2008 Smith et al. 2009a) which show 



that it is consistent with little or zero rotation |f] 

The resulting halo fractions are listed in Table [T] from 
which it can be seen that in general halo contamination 
is small. For all except the most metal-poor bins, the 
fraction of halo stars is less than 5 per cent and usually 
no more than 1 per cent. However, the lowest range in 
metallicity, which contains the largest fraction of halo 
stars, may be susceptible to halo contamination if the 
subtraction is imperfect. As expected, the number of 
halo stars increases as we move further from the plane. 
It should be noted that due to the spectroscopic selec- 
tion function of the SDSS survey, we do not expect these 
values to be representative of a volume-limited sample. 

3.2. Velocity Dispersions 

We use a maximum likelihood method to fit each vr 
and v z distribution as the sum of the disk and halo distri- 
butions, which are modelled using Gaussians. Although 
these velocity distributions are not expected to be ex- 
actly Gaussian, a maximum likelihood technique pro- 
vides more robust results than simply calculating the 
sample variance. The relative normalization is the halo 
fraction as calculated in Section |3.1[ and the velocity 
dispersions for the hal o are fixed at q-p = 138.2 kms -1 
and <t z = 89.3 kms -1 (Smith et al. 2009a I with no bulk 
motion. Note that in the calculation of the maximum 
likelihood, we take into account the errors on the indi- 
vidual velocities, as estimated in Section [2j 

We also use a maximum likelihood method to fit the v^ 
distribution. However, it is well known that the distribu- 
tion of ys for the disk is highly skewed and non-Gaussian 



'e.g.jg trombcrg 1927; Evans & Collett 1993] |Cuddeford 



i93| 

&; Binney||1994| |Binney & Merrineld||1998]). "Therefore. 



we require an asymmetric model to fit this distribution 
an d adopt the functional form described in equation (15) 
of Cuddeford & Binney ( 1994 ) , namely the distribution 



4 Our results have been found using the assumption that the halo 
is not rotating. However, for comparison purposes, we investigate 
the effect tha t a rotating halo wo uld have on our results (following 
the results of |Carollo et al.|2007| . This has a tendency to slightly 
increase our estimated halo tractions, but the effect on our results 
are small and our conclusions remain unchanged. Note also that 
the halo rotation is entirely degenerate with the assumed value for 
the local standar d of rest, which we have taken to b e 220 kms -1 
(see, for example, [Deason, Belokurov fc Evans|2011^ . 
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Fig. 1. — An example of one of the fi ts to our distribution of v^, 
using the method described in Section |3.2| This figure shows the 
most metal-poor bin for stars at heights of around f kpc below the 
plane (corresponding to the 10th column of Table [TJ . 



function for is given by, 



/(v )ocexp 



2al ( 



3 2/V4,/v c 



(1) 



x exp 



2oj 



Y <t> 



2v r 2 ln 



where y = 8 kpc/Rj. 

This results in a prediction for the distribution of v^ as 
a function of two free parameters (v c and Rd). To remove 
the dependence on vr in this equation we marginalise 
over the known distribution of vr for the disk stars in 
a given bin, as estimated using the method described in 
the previous paragraph. We use the observed v^ distribu- 
tion to calculate the likelihood distributions for v c and 
Rd- Given these likelihood distributions, we can then 
take the moments of equation ([I]) to estimate the mean 
and dispersion of v^ (and their uncertainties). Note that 
we do not infer any physical meaning from our determi- 
nation of v c and Rd - we are simply using the functional 
form given in equation ([I]) as a convenient description. 
We find that in most cases this provides a reasonable fit 
to the v^ distribution and should provide a more robust 
estimate than a sample variance. As before, we include 
a non-rotating Gaussian halo comp onent in the fit, with 
the dispersion fixed at 82.4 kms" 1 ( Smith et al.|2009a ). 
An example fit is shown in Figure 1. 



3.3. Covariance 

The covariance of the disk is calculated for each bin 
from the observed covariance using the relation 

(viVj)dh - fh(viVj) h 



(ViVj 



1-fh 



(2) 



where /h is the calculated value of the halo fraction, the 
subscripts 'dh' and 'h' refer to the combined disk plus 



halo sample and the halo respectively, and the subscripts 
i and j denote the coordinates R, 4> or z - In practice, we 
only analyze two cases, namely the tilt in the ellipsoid 
(i.e. (vrv z )) and the vertex deviation (i.e. (vrv^)). Note 
that when we calculate (viVj)dh we subtract the mean, so 
this should explicitly be written ((v; — (v;))(vj — (vj)))dh- 
For the halo stars, we assume that the ellipsoid is 
aligned in spherical polar coordinates, or in other words, 
it is pointing towards the Galactic centre. This has 
been shown to be a good approximation by |Smith et 



al. ( 2009b| ), who measured the offset to be no more than 
afew degrees for a sample of halo subdwarfs. 

Unlike the previous quantities, which were calculated 
using the maximum likelihood method, it is not possible 
to incorporate the errors on each individual velocity into 
this covariance calculation. To account for the effect of 
the errors, we instead adopt the following correction. For 
a particular star, the measured velocity, v; m , is equal to 
its true value, vj*, plus some displacement, Av;, which is 
due to the error in the measurement, i.e. Vj m = v^ + Av;. 
Given this, the true value of ViVj for this star is related 
to the measured value by 

(viVj)* = (v iVj ) m - Vi m Avj - Vj m Avi - AviAvj. (3) 

If we obtain the covariance by averaging this measure- 
ment over all stars in the bin, then clearly (vi m Avj) and 
(vj m Avi) should vanish. Therefore, we need only to cor- 
rect (viVj) by subtracting the term (AvjAvj), which we 
can approximate using our Monte Carlo estimates for the 
uncertainties (see Section [2]). 

This correction is important at high z, where the un- 
certainties on the tangential velocity are much higher. 
We find that this is especially problematic for the analy- 
sis of <tr z , since there are a larger fraction of stars at large 
R (and hence with large positive values of (AvrAv z )). 
In order to make our measurement more robust for <tr z , 
we reduce the radial range and exclude any stars outside 
of 7.5 kpc < R < 8.5 kpc. This reduces the number of 
stars in each bin (especially for large z, where some bins 
now have around 300 stars) and subsequently increases 
our statistical uncertainties, but we believe this is an ac- 
ceptable trade-off in order to reduce systematic errors. 

Once we have calculated the covariances, we use the 
following formula to obtain the angles 



#ij = 0.5 tan 



2<7ii 



(4) 



For both the covariances and the corresponding angles, 
we determine the uncertainties using a bootstrap tech- 
nique. We take 10,000 resamples from the observed dis- 
tribution of Vi and vj, with repetition, and calculate 
(v;Vj) for each resample, taking the dispersion of the re- 
sulting distribution as an estimate of the uncertainty in 
(v;Vj). The procedure is similar for fly, although in this 
case we also incorporate the uncertaint ies in the disper- 
sion, as calculated according to Section [3~2| 



4. KINEMATICS OF THE DISK 

All of the analysis in this section refers to properties of 
the disk population, recovered by subtraction of the halo 
contaminants. The values of the kinematic properties are 
summarized in Tabled] 



Slicing and dicing the Milky Way disk 



5 



TABLE 2 

Best-fit linear relations for the velocity dispersions. 



Component 


[Fe/H] 


cr(z = 0) 


da/dz 




(dex) 


( km s ) 


( If in o 1 lrl~if l^ 


^"R 


(— Fi ^fl 9*1 

^ — U.U, ^KJ'^iJ 


32.2 


11.3 


^R 


(_0.8, —0.5) 


35.3 


12.8 


CR 


(-1.5, -0.8) 


47.6 


7.9 




(-0.5, +0.2) 


22.4 


8.7 


<?<!> 


(-0.8, -0.5) 


25.9 


13.0 


<?<!> 


(-1.5, -0.8) 


32.1 


13.3 


&z 


(-0.5, +0.2) 


17.4 


9.3 


&z 


(-0.8, -0.5) 


28.9 


7.0 


0" z 


(-1.5, -0.8) 


38.2 


5.4 



4.1. Velocity Dispersions 

We plot the velocity dispersions as a function of z and 
metallicity in Fig. [2j The error bars represent the extent 
of the l-cr confidence interval determined by the maxi- 
mum likelihood fitting. We find that populations of lower 
metallicity are hotter than their more metal-rich counter- 
parts. This is a natural consequence of the correlation 
between age and metallicity - the metal-poor stars are 
typically older and hence have had more time to be af- 
fected by the heating mechanisms which cause disk stars 
to evolve away from the cold, circular orbits on which 
they were born. We also see that the velocity disper- 
sions increase on moving away from the Galactic plane, 
though there are hints that the velocity dispersion <r z 
may be saturating at high z. Table [2] reports the best-fit 
linear relations for these dispersions. Most are reason- 
ably well approximated by a linear fit, with the arguable 
exception of cr z which appears to saturate at large \z\. 

Such features have been found in a number of previ- 
ous studies. For example, Binney fc Merri fieldl (1998, 



sectio n 10.4.1) use a sample ot nearby stars ( |Stromgren| 



19871 to show how the velocity dispersions vary as a 
function of metallicity. These values are consistent with 
what we might expect if we extrapolate our results back 
to z = 0. More recent analyses have been done in the 
immediate solar neig hborhood using data from t he Hip- 
parcos satellite (e.g. Dehnen & Binney 1998a | or th e 
Geneva-Copenhagen survey (e.g. JN ordstrom et al.|2004| ), 
all of which are consistent with an extrapolation of our 
results. Analyses out of the plane have been limited 
due to the difficulty of obtaining reliable six-dimensional 
phase-space information for distant stars, but there are 
a num ber of such studies. For example, Spagna ct al. 
(2010) investigate trends in with both [Fe/HJ and z 



and Sbubiran et al. ( |2008| attempt to address the relation 
ilc 



between age and velocity dispersion for clump giants to 
wards the North Galactic Pole. The variation in er z with 
z, which is important as it can be used to trace the ver- 
tical potential of the disk (see Section [5]) , has rec eived 
a significant amount o f attention, most nota bly in |Kui- 
I jken fc Gilmore (1989) and more recently in Bond et al. 
( |2010| ). |Binneyj poiop has produced a distribution func- 
tion model which is able to predict the variation in <r z 
and we find that our results are in good agreement with 
this, although a detailed comparison is difficult due to 
the fact that we probe three separate metallicity ranges 
and do not provide an overall profile. 
Can we use the information contained in Fig. [2] to ad- 



dress the nature of the heating mechanisms? There are a 
variety of different phenomena that could act to heat the 
disk, from secular processes such as scattering due to spi- 
ral arms or mol ecular clouds (see section 8.4 of Binney & 
|Tremaine||2008| and references therein) , or external pro- 
cesse s such as accretion of sate llites onto the Milky Way 
(e.g. Villalobos & Helmi 2008). It has been postulated 
that the eccentricity structure of the t hick disk can be 
used to discriminate for mation models ( Sales et al.|2009 



Di Matteo et al. 2011). Recently a number of papers 
have i nvestigat e? this using ob servational data (Dierickx 
et al. 1120101 |W ilso n et al.||2011[ ). One of these works (|Di 
erickx et al.||2010) looks at the eccentricity distribution 



as a function of both height and metallicity and finds 
clear gradients, with the lower metallicity populations 
and those further from the plane are on more eccentric 
orbits; this matches the picture given by our data if we 
take Oft as a proxy for eccentricity 

These heating mechanisms may also result in differ- 
ences in the predicted anisotropics for the velocity ellip- 
soid. In order to investigate this, we plot the ratios of the 
dispersions in Fi g. [3| The solar neighborhood estimate 
for cr z /cr R » 0.5 |^len|T^f7| [Dehnen fc Binney||1998a| 
Nordstrom et al. 120041), is consistent with our metal-ricn 



sample. We detect a general increase in this ratio for 
lower metallicities, seemingly due to the fact that <r z has 
a stronger dependence on metallicity than or, i.e. a z in- 
creases faster than or as metallicity decreases, resulting 
in larger values for cr z /<7R. There is very little evidence 
for a gradient in z, and only a mild hint that the ratio is 
increasing with z. 

Theoretical studies of the expectations for cr z /<7R pre- 
dict a range of values less than 1. Classical work by 



Jenkins ( 1992 ) looked at heating by molecular clouds and 
spirals, and predicts ratios varying from 0.4 to 0.8 de- 
pending on the strength of the spirals, while a study into 
the effects of heating by molecular clouds alone predict s 



a value of around 0.6 (Ida, Kokuno & Makino 1983) 



Our results are in keeping with these predictions! In 
general, it is thought that scattering off spiral arms is 
more efficient at heating in the plane, while molecular 
clouds are more efficient at heating perpendicular to the 
plane. In this picture one might interpret the gradient 
with metallicity to claim that more metal-poor stars (i.e. 
older and hotter populations) are less affected by heating 
from spiral arms. Whatever the reason, it is clear that 
the metal-poor stars are significantly more isotropic in 
their kinematics than their metal-rich counterparts. 

Simulations of satellite accretion also predict similar 
ratios for t his ratio. For example, IVillalobos fc Helmi 
( |2008[ ) and IVillalobos, Kazantzidis fc Helmi| (|2010[ ) sin> 
ulate the formation ot the thick disk through heating via 
accretion, and also predict a wide range of values for 
Cz/cr (from ~ 0.4 to 0.9) depending on the parameters 
chosen to set-up their simulations (e.g. satellite mass 
ratio, orbital inclination, etc). Therefore, it seems that 
theoretical expectations cover a range of values, meaning 
that our data does not simply exclude a particular mech- 
anism. However, now that the velocity dispersions can 
be found as a function of height and metallicity, it should 
help to constrain future theoretical models by excluding 
specific ranges of parameter space. Very few models are 
able to predict the detailed behavior of these ratios and 
how they vary as a function of metallicity or z, which will 
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clearly become more important as models are refined. 

In recent years, there have been a number of attempts 
to measure <7 z /ctr for external galaxies. Even though 
such determinations are model dependent and based on 
a number of assumptions, it is interesting to compare 
these to the Milky Way. |Gerssen, Kuijken fc Mcrrifield 



( 1997 20001) looked at the large spiral galaxies NGC 488 



and NGC 2985, finding ratios of 0.7 and 0.8, respectively, 
i.e. slightly larger than the classical Milky Way dete rmi- 
nation of 0.5. Similar resu lts were also found by |van| 



der Kruit & de Grijs ( 1999 ) in their study of a compi 
lation of around 40 edge-on disk galaxies, with typical 
values ranging from 0.5 to 0.7 and with no evidence for 
a correlation between these values and the galaxy size or 
morphological classification (from Sb to Sd). Although 
larger than the Milky Way value, these ratios for external 
galaxies are consistent with our data; as can be seen in 
Fig. [3j as one departs from the plane and moves to lower 
metalhcities values of ct z /(tr, = 0.5 — 0.8 are typical. 

The ratio o$ / cxr also provides information about heat- 
ing processes, though this has received less attention in 
the literature. The value in the solar neighborhood is of 
interest as it is related to the slope of the rotation curve of 
the Galaxy. The ratio er^/cr^, which is sometimes refers 
to as Oort's ratio, is predicted to lie around 0.5 to 0.6 



depending on th e shape of the rotation curve ( Kuijken fe 
Tremaine||1991l see also Evans & Collett 1993). An ex- 



tre 



rapolation back to the solar neighborhood gives a good 
agreement with thi s expect ation, similar to the results of 
Dehnen & Binney ( 1998a[ ). In Fig. [3j we can see that 



the gradient with metallicity is less prominent than for 
Cz/cr, although it also looks like in general the metal- 
poor stars are more isotropic. As with <7 z /<7r, there is 
only weak indication that stars further from the plane 
exhibit larger ratios; in particular it is worth noting that 
the distribution for the metal- rich stars is essentially flat. 

4.2. The Rotation Lag 

It has long been understood that there is a correlation 
between the speed at which a population of stars rotates 
around the Galactic centre and the velocity dispersion of 
this population. This is known as the asymmetr ic drift 
(e.g. section 10.3.1 of Binney & Merrificld 1998) and is 



clearly evident in the solar neighborhood (e.g. Dehnen & 
Binney 1998a). The behavior out of the plane is more 
controversial. Whilst the existence of a correlation be- 



in numerous studies (e.g. 


Wysc & Gilmorc 1986 Majew- 


ski 


1992 


Chiba & Beers 


2000 Girard et al. 2000), the 



situation regarding trends with metallicity are less clear 
In recent years there have bee n at least two paper s with 
seemingly irreconcilable views; Spagna et al. ( |2010 ) claim 
to have detected a large gr adient, at a leyelot around 40 
to 50 kms^ 1 dex -1 , while Ivezic et al. (20081 claim an 



essentially flat relation with no detectable gradient. 

Our findings are presented in Fig. [4j where it is imme- 
diately evident that there are clear trends in the lag. The 
hotter populations (i.e. those with lower metallicity or 
those further from the plane) exhibit significantly more 
lag than their colder counterparts, with the metal-poor 
stars rotating more than 80 kms^ 1 slower than the LSR 
at 2 kpc from the plane. The gradient of the lag with re- 
spect to |z| varies from around 15 to 40 kms _1 kpc _1 , de- 
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Fig. 2. — Dispersions <tr, <t z and <T0 as functions of z and metal- 
licity. The triangles, squares and crosses correspond to metallic- 
ity ranges -1.5 < [Fe/H] < -0.8, -0.8 < [Fe/H] < -0.5 and 
-0.5 < [Fe/H] < 0.5, respectively. 

pending on metallicity. This is in agreement with many 
of the literature values which typically find values around 



Girard et al. 



2006). It is curious 



30 kms _1 kpc~ (e.g. 

to note that for the more metal-rich stars there appears 
to be a saturation in the level of the lag as one moves 
from the plane; from Fig. [4] it can be seen that while 
the metal-poor stars exhibit an almost linear trend, the 
gradient for the metal-rich and intermediate-metallicity 
populations becomes shallower for z > 1 kpc. 

Despite the strong theoretical grounds for expecting 
the lag to be larger for hotter (and hence more metal- 
poor) populations, as is seen in our data, some authors 
disagree with this claim. In particular, Ivezic et al. (2008) 
have argued that such a trend with metallicity is not ob- 
served in a sample of SDSS data. Th eir viewpoint is 
seemin gly lent weight by the findings of |Loebman et ah 
( |2010| , who analyzed an N-body+SPH disk galaxy sim- 
ulation and found no clear trend when plotting rotation 
velocity versus metallicity (see fig. 10). However, this 
does not imply that such a trend is entirely absent from 
their simulations; their fig. 9 shows that young stars do 
indeed possess larger rotation velocities compared to the 
older stars and als o that there is a clear age-metallicity 
gradient. If instead Loebman et al. ( 2010 ) slice their data 
in the same manner as we have (i.e. taken four bins in 
|z| from 0.5 to 2 kpc and three bins in metallicity), then 
they see behavior which is qualitatively identical to our 
Fig. |4]with a clear gradient in metallicity (V. Debattista, 
private communication) . 
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In the solar neighborhood, Dehnen & Binney (1998a I 
found that there is a linear relation between lag and 
the square of the radial velocity dispersion, which is in 
good agreement with t heory (e.g. section 4.8.2a of Bin- 
|ney fc Tremaine 2008). However, it is less clear what 
is expected as we depart from the plane. We show 
the observational result for our data in Fig. [5j As 
in the solar neighborhood, there appears to be a rela- 
tively tight correlation with crj| that is independent of 
metallicity. However, the gradient is much steeper than 
that of the solar neighborhood, where it is found that 
(v^) — v L sr = (v R 2 )/(80 ± 5 kms -1 ). The relation is 
reasonably well-fit by an quadratic relation, as shown in 
Fig. [5] where we find, 



(v*) - v LSR = 0.0149ct r 



1.21 x 10 -6 a R . 



(5) 



Note that our analysis h as been carried out using the 



Dehnen & Binney ( [l998a ) value for the Sun's motion 
with respect to the l ocal standard of rest. This has re- 



cently been revised by |Binney ( 2010 1 , with the velocity in 
the direction of rotation being increased from 5.2 kms -1 
to 11 kms -1 . If we adopt this newer value, then it im- 
plies that the lag which we measure will be too large by 
around ~ 5 kms -1 . However, since this offset should be 
applied uniformly to all of our sample, the gradients in 
lag with metallicity and height from the plane, which we 
can see from Fig. |4j will be unaffected. 

4.3. Radial and Vertical Bulk Motions 



Fig. 4. — The rotational lag, plotted against z. Velocities are 
rel ative to the Local Standard of Rest, assumed to be as determined 
by |Dehnen fc Binney (1998al, which means that in this plot the 
Sun would lie at (0, -5 kms -1 ). The triangles, squares and crosses 
correspond to metallicity ranges —1.5 < [Fe/H] < —0.8, —0.8 < 
[Fe/H] < -0.5 and -0.5 < [Fe/H] < 0.5, respectively. 



The bulk motion for the other two components, v z and 
vr are presented in Fig. [6] and Table [l] Curiously (v z ), 
whilst appearing consistent with zero for the metal-poor 
stars, seems to deviate from zero on moving away from 
the plane for the more metal-rich samples. Although 
tentative, this is an intriguing result and, if confirmed by 
further observations, could have important implications. 

In the lower panel of Fig. |6j we see that the radial com- 
ponent also seems to exhibit a net motion, with the sam- 
ple moving outwards from the Galactic centre . This has 
also been seen in data from the RAVE survey (Siebert et 
al.|2011 ), using a sample of stars with \z\ < 1 kpc. Mere, 
we seem to be detecting similar behavior out to \z\ ~ 2 
kpc, though there are no clear trends with either |z| or 



met allicit y. This is arguably surprising, since Siebert et 
al. ( |2011[ ) suggested that such behavior should be due 
to non-axisymmetric components in the Galaxy, such as 
the bar or spiral arms, which one might expect to have 
a greater effect on stars in the plane. The ot her poten- 
tial sourc e of non-axisymmetry postulated by Siebert ct 
al. poll) was that of ellipticity in the outer dark mat- 
ter halo, which we might expect to lead to a lack of a 
gradient with |z|. 

However, there are difficulties associated with measur- 
ing such quantities and various systematic biases could 
affect these results, for example mistaken assumptions 
regarding the Local Standard of Rest or the Sun's pecu- 
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Fig. 5. — Rotational lag as a function of radial velocity dis- 
persion. The triangles, squares and crosses correspond to metal- 
licity ranges -1.5 < [Fe/H] < -0.8, -0.8 < [Fe/H] < -0.5 
and —0.5 < [Fe/H] < 0.5, respectively. The dotte d line corre- 
sponds to the solar-neighborhood relation from[pchncn & Binney| 
(1998a |. The solid line denotes an empirical fit with the lag equal 
to 0.0149o-| + 1.21 X 10- 6 cr* . 

liar motion. 



4.4. Tilt of the Velocity Ellipsoid 

The covariance between v^ and v z , which is often re- 
ferred to as the tilt term, is an important quantity, which 
we can attempt to determine with our data. Little is 
known about how this behaves when one moves away 
from the Galactic plane. At z — 0, we expect this to van- 
ish due to symmetry arguments, but the behavior out of 
the plane is harder to measure as it requires distances to 
be determined to high accuracy. Despite these difhcul- 



understanding the shape of the potential (c.f. 


Smith et 


al.(2009b Binney & McMillan|2011 


Evans|20i: 


L}. Theo- 



between velocity ellipsoids parallel to the Galact ic plane 



to ones which point towards the Galactic c entre (Binney 



& Sperge ll|T M^|Kuiikerr^ Gilmore 1989| . The consen 
sus (e.g. Binney & Tremaine 2008]) seems to be that an 
alignment midway between cylindrical polar and spheri- 
cal polar is most reasonable. 

One of the recent determinations of the tilt was by 
Siebert et al. (20081. They used red clump giants from 
the HAVE survey to measure a tilt angle of 7°. 3 ± 1°.8 
at z ~ 1 kpc below the plane, with the orientation in- 
clined towards the Galactic centre. Their sample of stars 
covers a large range in z (from 0.5 to 1.5 kpc), which 



Fig. 6. — Mean velocities as a function of height from the plane. 
The triangles, squares and crosses correspond to disk stars with 
metallicities -1.5 < [Fe/H] < -0.8, -0.8 < [Fe/H] < -0.5 and 
[Fe/H] > —0.5, respectively. 

complicates the interpretation of their result, but it is 
clear that this value is close to what one would expect 
for an ellipsoid aligned with the centre of the Galaxy (i.e. 
tan _1 l/8 = 7°.l). This is an unexpected result as it lies 
at the extremum of the predicted values, but it appears to 
be robust as it has been supp orted by subsequent studies 
( Casetti-Dinescu et al.|[2011 ). 

Our results are given m Table [T] We have omitted 
the values for the metal-poor stars due to problems in 
accurately measuring their covariances. This is because 
it is very difficult to reliably correct for the halo con- 
tamination - for the metal-poor sample the covariance 
measurements are dominated by the contribution from 
the halo stars (since <7r z for the halo is much larger than 
for the disk) and so reliably extracting the disk covari- 
ance is practically impossible. We include the medium- 
metallicity sample in our analysis as the the halo con- 
tamination is smaller and hence our determination more 
robust. 

The covariances are plotted in Fig. [7] There appears 
to be a weak trend, with the magnitude of ctr z increasing 
slightly as we move away from the plane. If we now con- 
vert the covariances into an angles using equation (HI), 
we obtain the results shown in the lower panel of Fig. 
[7] The dotted line corresponds to what we would expect 
for a velocity ellipsoid aligned in spherical polar coordi- 
nates, and so we would expect our data to lie between 
this dotted line and the «r z = axis. Although there 
are large uncertainties, it appears that the metal-rich and 



Slicing and dicing the Milky Way disk 



9 





-100 
I? -200 

C\J 

E 

— -300 

N 

cc 
D 

-400 

-500 


-5 

CD 
CD 



cn 

CD 

N 

CC 



-10 



-15 



-20 



0.0 



0.5 



1.0 

|z| (kpc) 



1.5 



2.0 



Fig. 7. — The variation of ctr z and the corresponding angle or z , 
often referred to as the tilt, as a function of height from the plane. 
The dashed line is the assumed halo tilt (i.e. aligned in spheri- 
cal polars). The squares and crosses correspond to disk stars with 
metallicities -0.8 < [Fe/H] < -0.5 and [Fe/H] > -0.5, respec- 
tively. 

medium-metallicity stars are in general consi stent with 
the dotted line (and hence consistent with the |Siebert et 
al. 2008 result). A couple of points are more than one- 
sigma away from this line, but we believe this is probably 
an artifact; it is very difficult to explain such behavior 
for a disk in equilibrium, which implies that either our 
disk is not in equilibrium (due to accretion remnants af- 
fecting the distribution, or due to transient effects from 
resonances associated with spiral arms or the bar) or our 
data are suffering from observational biases. 

To conclude, although this is a difficult measurement 
with large uncertainties, it does appear that the tilt angle 
fo r the disk is surprisin gly large, reinforcing the findings 
of |Siebert etldT] p008| . 



4.5. The Vertex Deviation 

The orientation of the velocity ellipsoid in the (vr, v^,)- 
plane, often referred to as the vertex deviation, is also of 
interest. In the immediate solar neighb orhood, this has 
been found to be around 10° to 20° (|Dehnen & Bin- 

T00" 



ney 



km s 



1998a) 
7 



with corresponding covariance or^ 
For an axisymmetric system this term should 



vanish and the fact that it is non-zero is usually at- 
tributed to resonances from cither the Galactic bar or 



spiral arms 


(e.g. |Dehnen|| 1999| |Quillen||2003 


De Simone 


et al.||2004 


Minchev et al.||2U10f . 



ation out of the plane. Unlike the tilt angle discussed 
in the previous section, the vertex deviation is somewhat 
easier to measure. Since our data are located towards the 
South Galactic cap, velocities in the (vr, v^)-plane are 
mostly derived from the proper-motion data and, as a 
consequence, correlated uncertainties are less of a prob- 
lem. The halo correction is also easier to handle, al- 
though it still makes it difficult to robustly measure or^ 
for the metal-poor sample. Despite the fact that the halo 
= 0, the mean v^ is significantly offset from that of 
the disk and the dispersions are large. Therefore, unless 
the halo velocity distributions are well sampled (which 
they are not), there are likely to be large statistical fluc- 
tuations due to the influence of the halo contamination 
between different bins, even though the mean correction 
will be zero. 

Our results are presented in Fig. [8] and Table [T] As ex- 
pected, there are large fluctuations in the measurement 
of ctr0 for the most metal-poor stars For the rest, there is 
a clear p ositive signal, which is consi stent with the value 



found by Dehnen & Binney (|1998a 
on the most metal-rich stars 



If we concentrate 
which have much smaller 



errors than the more metal-poor stars, the distribution 
appears to be flat (assuming or^ = 100 km 2 s~ 2 at z=0) 
with a small increase for bin furthest from the plane. A 
similar rise is also seen for the medium-metallicity sam- 
ple. Instead of this weak positive gradient, we might 
naively expect that or^ should decrease with increasing 
height from the plane, since this offset is thought to be 
due to features which are stronger in the plane (namely 
the bar and spiral arms). However, this is only a very 
weak detection and so further study is required before 
any definitive statements can be made. 

5. AN APPLICATION: THE GALACTIC POTENTIAL 

There are many applications of the data we have pro- 
vided. We now show one such example where we use 
the vertical dispersions to compute the Galactic poten- 
tial close to the disk. 

For an isothermal disk it can be shown that, 



1 d{val) 
v dz 



<9$ 

~~dz~ 



(6) 



To obtain this relation we have neglected the contribu- 
tion of the tilt term (<7r z ); even if the tilt is pointing to- 
wards the Galactic centre (as suggested by Fig. [7j), then 
the contribution of these neglected terms to equation ([6| 
will be at a level of around five p er cent (see equation 
4.271 o f [Binney fc Tremaine|2008| ) . From equation ^ it 
is trivial to show that, 



v(z) = v(0) exp 



$(0) - $(z) 



(7) 



However, it is clear that the disk is not isothermal. If 
we generalize equation ^ to a population with disper- 
sion c z (z) then we obtain the following formula, 



Kz) = a z 2 (o) 

1/(0) a 



■ exp 



dz 



1 <9$ 

er 2 dz 



By taking a potential of the form, 



A$ = $(z) - $(0) = a\z\ + b\z\ 



(8) 
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Fig. 8. — The variation of <tj^ and the corresponding angle 
Q R0i often referred to as the vertex deviation, as a function of 
height from the plane. The triangles, squares and crosses cor- 
respond to disk stars with metallicitics —1.5 < [Fe/H] < —0.8, 
-0.8 < [Fe/H] < -0.5 and [Fe/H] > -0.5, respectively. 

and retaining only the first two terms, we obtain the 
following relation, 



v{z)o*(z) = K0)a z 2 (0) 



dzV(z') [a + 2bz'] . (10) 



This equation allows us to constrain the potential of 
the disk using our measured dispersions, provided we 
know the form of the density distribution for our tracer 
populations, i.e. our three metallicity ranges. Unfor- 
tunately this proviso is not an easy one to overcome. 
The complicated nature of the SDSS spectroscopic selec- 
tion function prohibits us from using our current spec- 
troscopic dataset to determine the density distribution. 
Therefore we resort to photometric number counts, us - 
ing the overall density distribution of Juric et al ( 2008 ) , 
convolved w i th the metallicity distribution funct ion of 
Ivezic et al. (2008 as revised in the appendix of Bond 
et al. pZUrD ). We cannot use these density distributions 
directly as they provide unrealistic predictions for the 
behavior at small z where, for the intermediate- and low- 
metallicity ranges, the resulting density actually rises as 
z increases (see the lower panel of Fig. [9]). To remedy 
this we assume that each of our t hree components follow 
a sech to the power 0.4 profile (Banerjee & Jog 20071 
and fit these functions to the Juric/ Ivezic profile over the 
range 0.5 < z/kpc < 2. For the range < z/kpc < 0.5, 
where the individual profiles are unreliable, we require 
the sum of our three density distributions to match the 



Juric total density distribution. We therefore include 
the scale-heights and normalizations of the three com- 
ponents as free parameters in our fit. Since the Juric 
disk profile consists of a thin- and thick-component, we 
similarly assume each of our metallicity populations con- 
sists of two components, keeping the scale-height ratio 
and the normalization of the thin- and thick-component 
fixed at the values adopted by Juric (zh,thick/zh,thin = 3 
and ^thick(0)Mhin(0) = 0.13). 

This means that we have a total of 11 free parameters: 
two to describe the potential plus three parameters per 
metallicity range (c z (0) and a normalization and scale- 
height for the thin-component of the density profile) . We 
then fit simultaneously the Juric/Ivezic den sity profiles 
and the dispersion profiles (via equation 10 1 by using a 
standard % 2 technique. We assume that the errors on 
the density profiles scale with y/v. 

Since our data do not constrain the dispersion profile 
for z < 0.5 kpc, we add an additional metal-rich data- 
point for the immediate solar-neighbo rhood using data 
from the Geneva-Copenhagen survey ([Nordstrom et al. 



2004| IHolmberg, Nordstrom fc Andersen|2009| . We take 
the 252 stars with [Fe/HJ > —0.5 dex, parallax error less 
than 13 per cent and distances less than 100 pc and find 
that the ct z = 15.1 ± 0.7 kms" 1 . 

As we are simultaneously fitting both the density and 
the dispersion profiles, it is in some sense arbitrary how 
we weight these two components of the fit. We have 
chosen to construct our fit so that the overall x 2 is ap- 
proximately equal for each of the two components, which 
implies that our fit is not dominated by either the density 
or dispersion profiles. 

The resulting best fit is shown in Fig. [9] with parame- 
ters given in Table [3j Although the Juric/Ivesic relations 
are not necessarily well represented by our sech profiles, 
note that the total density distribution from our three 
components provides a good match to the overall Juric 
et al (2008) distribution. However, it should be noted 
tfiat the fit relies strongly on our assumed Juric/Ivezic 
density profiles, which are subject to uncertainties that 
may bias our result. For example, one clear failing of our 
model is that the solar-neighborhood normalizations do 
not match observations. In particular ^(0) for the metal- 
rich sample is only around twice that of the intermediate- 
metallicity sample, whereas we know that in the solar- 
ncighborhoo d stars with [Fe/H] > —0.5 make up around 



95 per cent ( jHolmberg, Nordstrom fc Andersen 2009J) . 
As a consequence our resulting potential should not be 
over-interpreted. 

Given these caveats, it is reassuring to see that our 
potential (shown in the top panel of Fig. |9|) is in good 
agreement with existing models. Most notably our sim- 
ple ijnodej_^s_^n_e2ccOTtiojialagreement with the models 
of Dehnen & Binney (1998b hereafter DB), especially 
their Model 1. Although there is a slight discrepancy 
at small z (z < 300 pc), in general our potential matches 
these DB models to within 15 per cent out to 4 kpc. We 
have chosen to compare our potential to two DB models 

5 This discrepancy is most-likely due to the fact that in our 
simple model for the potential (equation^ one component of the 
mass distribution is an unphysical infinite razor-thin sheet. A bet- 
ter match at small z could be found if one distributed this mass in 
a more realistic manner. 
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(Model 1 and Model 4) because these can be considered 
to be two extreme cases of the models considered in their 
study, with Model 1 being the least halo-dominated and 
M odel 4 being the most halo- dominated (see section 2.7 



of |Binney fe Trem ainc 2008 for a detailed comparison 
of these two models, where Models I and II correspond 
to DB Models 1 and 4). Over the range where we have 
kinematic data (0.5 < z/kpc < 2) it appears that our po- 
tential favors Model 1, namely the model where the disk 
dominates the circular speed at the solar neighborhood. 

Once we have an estimate for A$ we can use this to 
probe the vertical mass distribution in the disk through 
Poisson's equation (9 2 $/9z 2 = AirGp). Our simple 
model corresponds to an infinite razor-thin sheet with 
a surface mass density of a/(27rG) = 32.5 M pc~ 2 , 
embedded in a uniform background with mass density 
6/(27rG) = 0.015 M Q pc~ 3 . Note that this background 
mass is close to the 0.014M Q pc -3 predicted usi ng isother- 
mal spherical h alo models (equation 4.279 of Binney & 
Tremaine|2008| . If we assume our background mass rep- 
resents the dark halo, it corresponds to a local dark mat- 
ter density of 0.57 GeV cm -3 , which is noticeably larger 
than the canonical value of 0.30 G eV cm~ 3 typically as 

alf 



sumed (e.g. 



various authors (e.g. 



Jungm an et al||199 6 
Gates ct al 



Ik 
P9l 



As pointed out by 



995 Weber & dc Boer 



2010 Garbari et al.|[2(jTT ), the local dark matter density 



is uncertain by a tactor of at least 2. Our analysis adds 
still more weight to the argument the local halo density 
may be substantially underestimated by the canonical 
value of 0.30 GeV cm -3 , and this is of immediate inter- 
est to dark matter experimentalists. 

Perhaps more robust than the local mass density is the 
surface mass density. By integrating our mass distribu- 
tion we obtain a total surface mass density of Si.ikpc = 
66 M0pc -2 , which agrees well with the classical value of 
71 ± 6 M pc- 2 from |Kuijken fc Gilmore| ( |l99l| ). If we 
integrate beyond 1.1 kpc, we find £ 2 k pc = 94 M pc~ 2 
and XUkpc = 155M Q pc~ 2 . 

Although we would ideally like to quote uncertainties 
for the above quantities, this is difficult due to the num- 
ber of assumptions and approximations in the adopted 
prescription. In particular, this method is dependent 
on what we assume for the density distributions for our 
tracer populations (in our case, the Juric/Ivezic profiles). 
As a consequence, the formal errors from the \ 2 fitting 
will be meaningless and we choose not to quote them. 
We have also experimented by adding a third term in 
the expansion of the potential (equation 9|, but found 
that gave no statistical improvement to the fit and had 
no significant influence on the resulting potential. 

6. CONCLUSIONS 

We have undertaken a kinematic study of the Galactic 
disk using data from the SDSS equatorial stripe region 
towards the South Galactic cap (Stripe 82) By combin- 



ing spectrosc opic data from Lee et al. 



( |2008[ ) with proper 

motions from |Bramich et al.' l 2008| ), we have constructed 
a sample of 7280 disk stars with full 3-dimensional posi- 
tions and velocities, along with [Fe/H]. This data have 
allowed us to probe the kinematics of the disk, tracing 
means and dispersions as a function of height from the 
plane and [Fe/H]. 
These data can be used to investigate the evolution of 




Fig. 9. — The results of the application to constrain the potential 
of the Galactic disk, as described in Section [5] The bottom panel 
shows the density distribution. In the range 0.5 < z/kpc < 2 
wc fit scch profiles to the metal-rich (cyan), mcdium-mctallicity 
(blue) and metal-poor (black) profiles derived from the work of 
Juric/Ivezic. In the range z/kpc < 0.5 we fit to the total density 
profile (thick black line). The total of our three density components 
is given by the thin red line. The middle panel shows the fit to 
the dispersion profiles, with the data taken from Fig. [2] plus an 
additional metal-rich point at z fa from the Geneva-Copenhagen 
survey. The upper panel shows the potential resulting from these 
fits. For the purposes of comparison, included in the upper panel 
are models for the potential taken from literature sources, namely 
Dehnen fc Binney| (1998 - Model 1, dotted; Model 4, dashed) and 
Fellhauer| (2UU6 - dot-dashed) . 



the disk, as stars which are born on circular orbits in 
the plane are heated by various mechanisms, such as the 
those caused by spiral arms, the bar, molecular clouds 
or accretion events. One avenue for addressing this is 
to investigate the ratio er z /<7R. For the metal- rich stars 
we measure this to be 0.6, which is consistent with pre- 
dictions for disk hea ting via spiral arms and molecular 
clouds (|Jenkins|1992[ ). However, predictions from models 
of satellite accretion also cover the range we observe and 
so it is difficult to deduce any strong conclusions. We also 
investigate the ratio o^/ctr, and find that in general stars 
further from the plane exhibit larger ratios, although this 
trend is less evident for the more metal-rich stars. It will 
be interesting in future for models to attempt to explain 
the gradients, both in z and [Fe/H], which are evident in 
our data. In particular, it is clear that the metal-poor 
stars are more isotropic in their kinematics than their 
metal-rich counterparts. 

Our data also allow us to measure the covariances and 
constrain the orientation of the velocity ellipsoid. We 
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TABLE 3 

Best-fit parameters from the potential fitting procedure. 
Parameters a and b correspond to the potential given in 

EQUATION |9]|, 1/(0) IS THE SOLAR-NEIGHBORHOOD DENSITY 
NORMALIZATION AND Zh,thin AND Zh. thick ARE THE SCALE-HEIGHTS 
OF THE THIN AND THICK COMPONENTS OF THE DISK PROFILE, 
RESPECTIVELY. NOTE THAT THE SCALE-HEIGHT OF THE THICK 
COMPONENT IS NOT A FREE PARAMETER AND IS KEPT FIX ED AT 
THREE TIMES THAT OF THE THIN COMPONENT (FOLLOWING |JURIC~| 
I ETAL|2008) . 



1 til CLI111_- Lt_.l 




/ 

R ti vi cm 




LI {r^lLl rvL/L, 1 






-876.7 


b (km 2 s -2 kpc -2 ) 






-414.4 


i/(0) (10- 3 M Q pc- 3 ) 


(" 


0.5, +0.2) 


6.1 


Zh,thin (kpc) 






0.05 


Zh, thick (kpc) 






0.16 


<j R (0) (kms- 1 ) 






14.9 


i/(0) (10- 3 M Q pc- 3 ) 


(" 


0.8,-0.5) 


2.7 


Zh.thin (kpc) 
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0.14 


Zh, thick (kpc) 






0.42 


<7 R (0) (kms -1 ) 






26.3 


1/(0) (10- 3 M Q pc- 3 ) 


(" 


1.5,-0.8) 


0.8 


Zh.thin (kpc) 






0.26 


Zh, thick (kpc) 






0.78 


<t r (0) (kms -1 ) 






38.2 



found that the tilt term (q;r z ) is not consistent with zero, 
i.e. the ellipsoid is not aligned with cylindrical polar co- 
ordinates. Our results are consistent with previous stud- 
ies, which have found that the vertical component of the 
ellipsoid is clo se to being aligned in spherical polar co- 
ordinates (e.g. Siebert et al.|2008 ). The vertex deviation 
(oiH(f>) is found to be consistent with an extrapolation 
to the solar-neighborhood, with a marginal detection of 
an increase as one moves away from the plane. If this 
gradient can be confirmed it may have interesting impli- 
cations, as the mechanisms that are believed to drive this 
term away from zero (e.g. the bar or spiral arms) should 
have greater influence in the plane. 

In order to address the nature of disk heating and to 
disentangle the contributions from various mechanisms, 
one needs to go beyond the work presented here. The 
most crucial improvement will be the ability to make 
direct estimates for stellar ages, rather than relying on 
correlations with metallicity. Ages are notoriously dif- 
ficult to measure robustly, even for surveys of bright 



nearb y stars (e.g. Holmbcrg, Nordstrom & Andersen 
20091. Extending such studies beyond our immediate 
solar-neighborhood will be a difficul t task, but one tha t 
is currently showing promise (e.g. Burnett et al. 2011). 



One aspect that will help us in this effort is by told- 
ing in measurements of alpha-element abundances that 
are now b eing determined routinely for vast numbers of 



stars (e.g.|Bovy et al.||2011[ |Lee et al.pOlT 



al. 2011 Ruchti et al. 2010). In particular Bovy et al. 
( 2T)lip ana Lee et al.| ( |201l[ ) show what is possible with 



Navarro et 



SUSS data when one incorporates information on alpha- 
element abundances. 

In recent years there has been some debate regarding 
the rotation lag of the disk. The lag has been known 
about for many years, wit h recent determin ations of 
around 30 kms _1 kpc _1 (e.g. Girard et al. 20061. We find 



gradients of around 15 to 40 kms _1 kpc _1 . These gradi- 
ents depend on height from the plane (with the gradients 
in general becoming shallower as z increases) and also on 
metallicity. Our findings address one particular bone of 
contention, which is whether there exists a correlation 
between rotation lag an d metallicity (e.g. |Ivezic et al. 
|2008[[Spagna et ah 2010 ). Our data clearly show that the 
more metal-poor (and hence hotter) stars exhibit greater 
lag than their metal-rich counterparts, with the metal- 
poor stars reaching a lag of more than 80 kms -1 at 2 
kpc from the plane. Our findings are in good qualitative 
agreement with simulations, although further work needs 
to be done to fully exploit these observational results. 

In passing, it is interesting to compare our results on 
the vertical gradient of the lag to that for neutral HI 
gas. For the Mil ky Way this has been found to be 15 ± 
4 kms _1 kpc~ 1 (Marasco & Fraternali 20111, which is 
comparable to our findings for metal-rich stars. This 
measurement is arguably easier to determine for external 



30 kms-^pc" 1 ( 


Heald et al. 2007 |Kamphuis et al. 2007 


Zschaechner et al 


. 201 1\. Although the gas response to 



the collisionless stars, the magnitude of the lags appear 
similar. 

Another limitation of our analysis is that we are con- 
fined to the 250 square degrees of Stripe 82. Although we 
are able to probe gradients with height from the plane, 
we are unable to draw any conclusions regarding radial 
gradients. Furthermore, the limited sky coverage means 
that we are susceptible to bias should kinematic sub- 
structure be present within any of our bins. Although 
no such substructure is immediately evident in our data, 
we cannot exclude such an occurrence. By extending 
this to the whole SDSS f ootprint using the r ecently re- 



leased 8th Data Release (Aihara et al. 20111 it will be 



possible to overcome some of these limitations, although 
the lack of radial coverage will persist as SDSS has very 
few fields at low Galactic latitude. Future work that will 
help to fill in the missing information at low latitudes 
includes the LAMOST telescope (e.g. |Wu et al.|[MT | 
surveys, which will in-part focus on the Galactic disk. 
Looking further ahead the field will undoubtedly be rev- 
olutionized by the Gaia satellite mission, which will pro- 
vide precise distance s and proper motions for billions of 
stars in our Galaxy ( Perryman et al.||2001 Sumi|[2009 ). 

We have also shown an application of Dur data, em- 
ploying the Jeans equations to provide a simple model 
of the Galactic potential close to the disk. We use the 
density model of |Juric et al| ( |2008[ ) and our measured val- 
ues on the variation of cr z . Our model does a good job 
of reproducing the kinematic data for stars in all three 
ranges of metallicity. T he model is in excellent ag ree 



ment with Model 1 from Dehnen & Binney (1998b I, in- 
dicating a preference for a model where the disk, rather 
than the halo, dominates the circular speed at the so- 
lar neighborhood. The main uncertainty in this study 
is the underlying profile of the tracer populations, since 
this is something we cannot recover from our data due 
to the complicated selection function of the SDSS spec- 
troscopic survey. Future unbiased spectroscopic surveys, 
such as those undertaken by the LAMOST telescope 
or the Gaia mission, will allow much more robust con- 
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straints to be made. A knowledge of the Galactic po- 
tential near the disk is clearly great importance. In par- 
ticular, it allows us to constrain the local dark matter 
density, which is a crucial piece of knowledge when eval- 
uating the pr ospects for direct and indirect d etection of 
dark matter ( |Bertone, Hooper fe; Silk||2005| and refer- 
ences th erein). Despite slow progre ss since the classical 
work of Kuijken & Gilmore (1991), this field is set to 
have a great impact in the coming years. 
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APPENDIX 

THE TURN-OFF CORRECTION FOR THE PHOTOMETRIC PARALLAX RELATION 



Our distances are estimated using the photometric parallax relation of Ivezic et al. ( 2008 1 , with one minor modifi- 
cation regarding the treatment of the main-sequence turn-off. Ivezic et ST. ~{ 2008 ) based their turn-off correction on 
the color-magnitude sequence for cluster M13, which has an age of around 10 Gyr and [Fe/H] = —1.54. Although this 
is a sound approach for stars belonging to the Galactic halo, disk stars will be (in general) both younger and more 
metal-rich than this cluster. Therefore we adopt an altern ative approach to modelling the turn-off correction. 

We follow the approach described in Smith et al.| ( |2009a[ ) , where stellar models were used to construct a metallicity- 
dependent turn-off correction. Since we cannot discriminate between different populations on a stars-by-star basis, we 
must construct a global correction which reflects the relative numbers of thin-disk, thick-disk and halo stars at a given 
[Fe/H]. This is done by constructing a toy model to represent the properties of our sample, with parameters given in 
TableH " 



The final ingredient for this calculation are theoretical stellar models, which we take from |Dotter et al. (2008). We 
downloaded t hree sets of models, c orresponding to each of our three Galactic components. The prescription we adopt 
is the same as Smith et al. (2009a). In brief we we shift each sequence so that it has M r = for (g—i) = 0.6 and then, 
for a given [Fe/HJ, we calculate a weighted mean M r as a function of (g—i) in the range 0.3 < (g—i) < 0.6, considering 
only model data up to the main-sequence turn-off. The weights are determined from the above toy model, so that we 
account for the contribution of each Galactic component at a given metallicity. 

We the n calculate the offset between the weighted mean model magnitude and the uncorrected relation of |Ivezic"et] 
al. (2008), i.e. given in equations (Al— A5). This is shown in Fig. 10 Clearly, for met allicities below sol ar all of the 



models lie around AM™ = 0, i.e. they agree with the uncorrected parallax relation of jlvezic et al.| ( |2008[ ). Therefore 
we decide not to incorporate any turn-off correction when calculating our distances. 

The only exception from this behavior is the solar-metallicity model, which is systematically offset by between 0.1 
and 0.2 magnitudes. We believe this is due to problems with the normalization of the models (which we do by forcing 
all models to pass through the point M r = at (g — i) = 0.6). In any case, this should have little effect on our 
results as very few metal-rich stars are blue enough to lie in this turn-off region; only one per cent of our sample have 
[Fe/H] > -0.25 and (g-i) < 0.6. 

There will be additional uncertainties on the parallax relation in this turn-off regime due to scatter in this correction. 
To account for this we calculate the standard deviation in M r when we calculate the mean. The scatter varies as a 
function of color and metallicity, but if we take the relation with the largest scatter (corresponding to [Fe/H] = — 1 
dex) we find the following relation, 

S (AM™) = 0.272 — 0.454 (g—i). (Al) 

When estimating distances for our stars we add this uncertainty in quadrature to the other sources of error for stars 
in this color range. 
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TABLE 4 

Parameters of the toy model used to construct the turn-off correction to the photo metric paral l ax re lation. The 

PARAMETERS ARE CHOSEN TO BE CONSISTENT WITH THE FOLLOWING REFER ENCES: DENSITIES FROM JURIC ET Al| 1 2008| , WHICH WE 
EVALUA TE AT Z = 1 KPC, I.E. THE M EAN VALUE OF OUR S AMPLE; [Fe/ HJ FROM jS oi'BIR AN ET AL. (2003| FOR THE DISK AND |f VEZIC ET AlT| 



(2008) FOR. THE HALO; [a/Fe] FRQm |VeNN ET AL~] | |2004| ; AGES FRO M |NORDSTROM ET AL.| | |^UL)4 FOR THE DISK AND CHOSEN TO MATCH 

THE HALO AGE USED IN ( |SMITH ET AL.|2U0W| 
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TABLE 1 

Kinematic properties of the disc. The values of <tj^ z and q^ b for stars with [Fe/H] < —0.8 dex are unreliable due to the difficulty in making the halo correction and hence have been omitted 

FROM THIS TABLE. NOTE THAT erj^ z AND aj^ z WERE CALCULATED USING ONLY STARS WITHIN 7.5 KPC < R < 8.5 KPC AND SO THE NUMBER USED IS LESS THAN THAT QUOTED IN THE SECOND COLUMN, ESPECIALLY AT 

LARGE Z. 



[Fe/H] (z) No. Halo <v R > 

(dex) (kpc| stars Frac. (kms - 1 ) 



( kms -1 ) 



< v z> "R <*4> "* »E. »Rz 

(kms -1 ) (kms -1 ) (kms -1 ) (kms -1 ) (km 2 .- 2 ] (°) 



"R<t> 
(km 2 s- 2 ) 



<*R</> 
(°) 



-0.5, +0.2) 
-0.5, +0.2) 
-0.5, +0.2) 
-0.5, +0.2) 

-0.8, -0.5) 
-0.8, -0.5) 
-0.8, -0.5) 
-0.8, -0.5) 

-1.5, -0.8) 
-1.5, -0.8) 
-1.5, -0.8) 
-1.5, -0.8) 



-0.59 736 

-0-86 736 0.00 

-1.11 736 0.00 

-1.53 736 0.02 



-0-59 543 
-0-96 543 



-1.26 
-1.69 



-1.32 
-1.73 



543 
543 



-0.60 541 
-0-99 541 



541 
541 



0.00 -O.O^H 
— 1.6 

+ 1.9 



0.00 
0.00 
0.01 
0.04 

0.11 
0.14 
0.17 
0.27 



3.9^ 
5.6^ 



„ + 2.4 
8 -2.4 



0.5 
6.1 
4.3 
4.0 



+ 2.0 
2.0 



+ 2.4 
2.3 



+ 2.6 
2.5 



+ 3.2 
3.2 



.+3.0 
-2.9 



4.3 
5.1 
7.6 



+ 3.2 
3.3 



+ 3.9 
4.0 



+4.6 
4.7 



9 +l.l 
-1.1 



-188.4 
-183.9 j 
-179. 3 4 



+ 1.2 
1.2 



g+2.1 

-2.1 



+2.1 
J -2.1 



6+2.6 
°-2.6 



4+3.3 
-3.3 



B +4.0 
°-4.0 



7 +0.9 
-1.0 



-1.9 
-6.0 



+ 1.1 
1.1 



+ 1.2 
1.3 



a+1.5 
8 -1.4 



+1.5 
" — 1.5 



„ + 1.8 
8 -1.8 



! + l-9 
1 -1.9 



0.6 
0.1 
1.0 
1.2 



+ 2.2 
2.1 



+ 2.6 
2.6 



+ 2.8 
2.8 



+ 3.4 
3-3 



38.8 
42.6 
43.5 
50.2 



+ 1.2 
1.3 



+ 1.5 

1.4 



+ 2.1 

2.2 



7+1.5 
-1.6 



9 + 1.9 
2 -1.8 



,+2.1 
-2.1 



6 +2.7 
-2.7 



52.9 
53.0 
61.3 
60.2 



+ 2.7 
2.7 



+ 3.2 
3.2 



27.7 
29.8 
33.6 

35 



+ 0.4 
0.4 



+ 0.1 
0.1 



3 + 1.0 
3 -0.9 



32.8+^0 



6 + 1.9 
-1.9 



2 + 2.4 
-2.5 



40.8 
42.6 
50.7 
55.3 



+ 2.0 
2.0 



+ 2.6 
2.7 



+ 2.7 
2.7 



+ 3.0 
3.0 



22.4 
25.7 
28.5 
30.7 



+ 0.7 
0.7 



+ 0.9 
0.9 



+ 1.0 
1.0 



+ 1.2 
1.2 



„ + 1.4 
8 -1.4 



, + 1.5 
-1.5 



+ 1.8 
-1.8 



40.2 
46.1 
45.7 
45.4 



+ 1.8 
1.8 



+ 3.0 




_+36-8 
-36.8 



-242.3+?^'? 



+ 104.2 
J -106.0 



-193.6 



+ 71.9 
70.9 



-154.0 
-181.0 
-311.3 



+ 127.6 
127.8 



+ 133.7 
135.0 



+ 179.3 
179.8 



-5 6 + 2 ' 1 
— 2.1 

-14.6± 2 -» 
-12. 0" 1 " 



-12.4 



+ 3.5 
3.5 



-12 7+ 4 ' 6 
12 - -4.6 

- 1 O 2+ 8,3 



110.4 
206.9 



+ 45.6 
44.7 



+ 60.9 
61-5 



100.3 



+ 76.8 
76.5 



+ 126.5 
-127.1 



Q + 144.6 
-144.0 



„ + 227.6 
-225.0 



B + 418.9 
-406.6 



„ + 473.9 
-471.1 



„ + 555.2 
-551.9 



, + 496.2 
-492.8 



8 3 + 34 
83 -3.4 

1 9 O 4 " 3- ^ 

7 4+ 6 - 5 
' -5.6 

16 a+ 4 6 



10-0 



+ 5.6 

5.5 



-6.6 
13.5 
18.3 



+8.3 
8.3 



+ 13.5 
13.5 



+ 7.3 
7.1 



, + 17.7 
-18.7 



„ + 18.9 
°-21.0 



-+36.8 
-32.1 



